Self-energy of a nodal fermion in a d-wave superconductor. 
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We re-consider the self-energy of a nodal (Dirac) fermion in a 2D d— wave superconductor. A 
conventional belief is that Im £(w,T) ~ max(Lu 3 , T 3 ). We show that k, T) for k along the 
nodal direction is actually a complex function of uj,P, and the deviation from the mass shell. In 
particular, the second-order self-energy diverges at a finite T when either uj or k — kp vanish. We 
show that the full summation of infinite diagrammatic series recovers a finite result for E, but the 
full ARPES spectral function is non-monotonic and has a kink whose location compared to the mass 
shell differs qualitatively for spin-and charge-mediated interactions. 



The physics of nodal quasiparticles in a 2D d— wave 
superconductor attracted a considerable attention since 
the early days of high T c era Q, 0, because of its uni- 
versality Q| and a relation to field theoretical studies of 
Dirac ferrnions . In the normal state of the cuprates, 
ARPES and other studies have found that the imaginary 
part of the fermionic self-energy along the nodal direc- 
tion, JmS(o;,T) 1 is roughly linear in both frequency and 
temperature 0, S In the superconducting state, a 
generic belief is that im£(cJ,T) should become smaller 
as whatever the mechanism is for the quasiparticle scat- 
tering in the normal state, it must be weakened in the 
superconducting state due to the gap opening Q ■ Recent 
ARPES results do show indeed that once bilayer splitting 
is resolved, the measured Jm£(a;,T) decreases below T c 
faster than linear 0, 0, EJ • 

In an s— wave superconductor, ImY>{uj,T) at uj < 3A 
vanishes at T — and is exponentially small at finite 
T In a d x 2_ y 2— superconductor, the gap vanishes 

along the diagonal directions in the Brillouin zone. Hence 
the scattering into low-energy states near the nodes gives 
rise to power-law uj and T dependences of ImY>{uj, T) at 
the lowest u> and T. 

Several groups have previously computed 0, E3, E-l 
the quasiparticle lifetime of a nodal fermion in a BCS 
d— wave superconductor at T = and at k = kp, and 
found the cubic frequency dependence 7mE(w,T = 0) oc 
ui 3 . A common belief was that (i) this dependence sur- 
vives at arbitrary ratio of uj and = vp(k — kp) along 
the nodal direction, and (ii) for a clean system, the tem- 
perature dependence is the same as the frequency depen- 
dence, i.e., 7mS(w,T) scales as uj 3 or T 3 , wichever is 
larger. 

In the present communication, we dispute this common 
belief. We show that the perturbative self-energy for a 
nodal fermion actually scales as T 7 ^ 2 / yfuj at T » uj and 
generic uj/tu and logarithmically diverges at the mass 
shell, u — £fc. This singular behavior has been largely 
missed in earlier studies (with the exception of , where 
T 7 / 2 / v /7 J has been reported). We show that the origin 
of this singular behavior is the same as of the mass-shell 



singularity in a 2D Landau Fermi liquid (LFL) [HDS- 
L/ike in a LFL, the singular behavior of the self-energy for 
a Dirac fermion is eliminated once infinite series of diver- 
gent terms are summed up, but the resulting spectral 
function has a kink at some finite uj — We found that 
the kink is located at > uj > efc, or at > ej, > uj de- 
pending on whether the effective interaction is in the spin 
or in the charge channel. From this perspective, detailed 
ARPES studies of nodal spectral function can qualita- 
tively distinguish between the theories for the cuprates 
based on spin or charge fluctuations. 

We consider a model of a d— wave superconductor with 
a quasiparticle dispersion 
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|k|, where 

k = (upfcx) va^ii), and k± and fey are deviations from 
the nodal point transverse to and along the Fermi sur- 
face , respectively (e^ = vpk±). The fcrmion-fermion 
interaction has both charge and spin components, U c (q) 
and U s (q). We explicitly verified in lengthy calculations 
that only intra- nodal interactions, governed by U c (0) and 
U s (0) contribute to the self-energy along the nodal di- 
rection, and we consider only these terms. Inter-nodal 
interactions do contribute to the self-energy away from 
the nodal direction, which we did not consider. 

We first re-analyzed the self-energy to the second order 
in the perturbation. We carried out calculations both in 
the Nambu formalism and using normal and anomalous 
Green's functions, and obtained identical results in both 
approaches. In the Nambu formalism, the mean field 
action near a given pair or nodes at kj? and — ki? is de- 
scribed by the four-component massless Dirac ferrnions 
(Nambu spinors) and has the following Lorentz invariant 
form: 



C = iV^Todr + v F -fid x + v A j 2 d y )^ 



(1) 



where if,+ (k) = (V>J(*), ^(~k),i>{ (A), Vi(-*)J , i> = 
ip + lo, 7o = I ® <J y , 71 = 1 cr x , 72 = a z <g> a z and a a are 
the Pauli matrices. The free-fermion Green's function is 
G = (ipip) = ^IfiP^/p 2 - The fermion- fermion interaction 
is the sum of the bilinear products of charge density and 
spin density operators tp^ipa and ijj^Gapijjp, respectively. 
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The perturbation theory is constructed in the same way 
as in LFL. 

In the conventional formalism, which we will follow 
below, the self-energy is obtained by summing up contri- 
butions from normal and anomalous Green's functions of 
intermediate fermions. The expression for ImT,(uj, k).T) 
obtained in this formalism is presented in We ver- 

ified and confirmed their result. The authors of Q ana- 
lyzed the self-energy numerically and argued that it fol- 
lows max(uj 3 , T 3 ) behavior. We evaluated the self-energy 
analytically and found new, singular behavior at uj ~ 
at T = 0, and at w,ej;<Tata finite T. To understand 
the origin of the singularity and to establish connection 
to the earlier work on the mass-shell singularity in LFL, 
it is convenient to view the self-energy as a convolution 
of the fcrmionic and particle-hole propagators: 

9 f 9 / ^ UJ + fl\ 

ImY,(uj, k, T)<xU 2 d 2 qdn ( coth — - tanh ^ j 

Imxo(£l,q)ImGo(w + fl,k + q) (2) 

where Go is the BCS fcrmionic Green function. In non- 
Nambu notations, the polarization bubble is the sum of 
normal and anomalous components, and is given by 



xo(n,?)« / d z i 



.E. 



I 2 + E+E_ 



E+E- (E, 



(3) 

where E± = |l±q/2|. At small << q, Imxo(fl,q) oc SI, 
as has already been established before |l2j. The convo- 
lution of this form with the linear frequency dependence 
of the fermionic density of states in a d— wave supercon- 
ductor gives rise to JmE(w, ek, T) oc (uj 3 ) at uj ~ T and 
generic e%/uj, in agreement with [7J. The behavior of the 
self-energy near the mass shell is, however, determined 
by |fi| w q, where Xo(^, (?) is singular. The singularity 
originates from the integration in © over I < q/2, and 
over directions of 1 which are nearly parallel to q. In this 
range, E+ + E- ~ q, i.e., the denominator almost van- 
ishes in a finite range of internal I, while the numerator 
in JSJ remains finite. Expanding E + + E- in the angle 9 
between 1 and q as E+ + E_ m q[l + l 2 9 2 /(q 2 - 4/ 2 ) + ...], 
substituting into |PJ and integrating over 0, we obtain 



Xo(n,g) =A 



vV - ($i + isy 



(4) 



where A > 0. We see that Imxo(fl,q) is nonzero when 

|fi| > q, and it diverges at — ► q. 

Substituting Imxo(Q,q) and IitiGq(uj + 0, k + q) oc 
\uj + fi| 5((ui + tt) 2 - (k + q) 2 ) into ©, we find that on 
the mass shell (i.e., when at uj = k), the position of the 
branch cut singularity in 7mx(f2, q) and the location of 
the 8— function in ImG(uj + Q, k + q) coincide provided k 
and q arc cither parallel (for uj' > 0) or antiparallel (when 
lu' < 0). Integrating near these directions, we obtain the 
singular part of the self-energy. At T = 0, we find that 



/mS(u, Ck,T — 0) is non-zero only when uj < ek (in 
conventional variables where Ek = |ejt|, and for uj < 0, 
relevant to ARPES experiments), and is discontinuous 
at uj = ek (Re's diverges logarithmically at this point). 
Explicitly, 



IiriE(uJ « ek, T = 0) = 
2tt 3 , 9 „ f vf \ 2 w 3 



E 2 



'(-'-M) (5) 



where Ef = VFkp/2 and u c , s = rnU c , s (0) / (2tt) 
At a finite T, we obtained in three limits: 



JmS(w = 0,k,T) 
ImE(w,k = 0,T) 



n7/2 



Te 2 -\e k \/2T 



r,7/2 



M' 



3/2 



M «r, 

M »t. 
M « r, 

M >> r. 



ME(w « e fe , T) ~ M 3 0(e fe - w) + 
' = (aT 2 + few 2 ) log | — | (6) 



e k -oj 



where a, b = 0(1). We see that at finite T, the second- 
order self-energy diverges logarithmically on the mass 
shell, and also diverges as 1/ \/\oj\ when both uj and ek 
tend to zero 0]. This result is very different from a 
simple T 3 behavior. The latter is only recovered when 
|oj| ~ T, and the system is at some distance away from 
the mass shell. 

As we have said, the singularity in the self-energy 
for a Dirac fermion is similar to the mass-shell singu- 
larity in the fermionic self-energy a 2D LFL [lH Il7| . 
In both cases, the singularity originates from the fact 
that the the pole in the fermionic Green's function al- 
most coincides with the branch-cut in the polarization 
bubble |lt|. There exists, however, an important dis- 
tinction between the two cases, which makes the singu- 
larity for Dirac fermions stronger than in a Fermi liq- 
uid. In a liquid, the quasiparticle dispersion changes 
sign at the Fermi surface, i.e., ek F + q = qcos9 can be 
both positive and negative. Accordingly, ImT,(uj,T) oc 
fd 2 qd6doj' q S(oj' + (uj - k) - qcosO)/^^') 2 - q 2 . and 
the angular integration reduces the 5 function to another 
square-root, so that the subsequent momentum integral 
only logarithmically depends on the uj — k, i.e., on the 
distance from the mass shell. For Dirac fermions, the 
dispersion Ek = \k\ is positive, and Ek F + q oc \q\ > 0. 
The angle integral still accounts for the logarithmic sin- 
gularity at uj = Ek, but in addition, the combination 
of the 5 function from the fermionic propagator and the 
square-root singularity in the polarization operator gives 
rise to the extra (T/Iuj]) 1 ^ 2 factor in the self-energy. 



3 



The discontinuity of the self-energy on the mass shell 
at T = and the divergence at T > imply that higher- 
order diagrams for the self-energy may also be relevant. 
Evaluating higher-order self-energy diagrams with extra 
particle- hole bubbles, we find that they actually diverge 
at the mass shell already at T — 0, and the divergences 
proliferate with the order of perturbation. Like in a 
LFL, the most divergent diagrams form ladder (RPA) 
series 0| > such that the full self-energy can be still repre- 
sented by Eq. (j2J), but now Xo(&>, <?) has to be replaced by 
the full susceptibility x(fi, g). This computational pro- 
cedure is justified when u c>s <C 1. 

At this stage, it becomes essential whether the effective 
interaction is in the charge or in the spin channel, as the 
RPA renormalizations in charge and spin channels have 
opposite signs. Summing up the bubble and the ladder 
diagrams, we obtain the following expressions for the full 
Xc(ft,q) and x a (n,g) 



/mE(w, k) becomes continuous. However, it still vanishes 
at > to > e k . We found 



Xo{M,q) 



i + u cXo (n, q ) 



Xs(ci,q) = 



Xo{M,q) 



l-E7.xo(fi,?) 



(7) 

where xo(^><?) is given by Q (recall that C/ C , s > and 
Re X o(Q,q) > 0). For a charge-mediated interaction, 
there is no pole outside the particle-hole continuum. Ac- 
cordingly, the RPA renormalization only eliminates the 
square-root singularity in Im X (£l, q) at Q = q. The full 
Imx(£l, q) then behaves as in Fig. [IJi: it initially fol- 
lows Imxo(Q, q) and increases at approaching Q = q, but 
passes through a maximum and vanishes at = q. For 
X s (f2, q), the sign of the renormalization is different, and 
the interaction not only eliminates the square-root diver- 
gence at 51 = q, but also generates a zero-sound-type pole 
outside the particle- hole continuum, at £1 = q(l — Bu 2 q 2 ), 
B ~ l//c| > 0, as on Fig. Hh- 



FIG. 1: The behavior of the imaginary part of the full bosonic 
propagator Im x(^i <l) f° r interacting Dirac fermions. Panels a 
and b are for the interaction in the charge and in spin channel, 
respectively. For spin interaction, x(^i <?) develops a zero- 
sound pole outside the particle-hole continuum. The dashed 
line is the imaginary part of the bosonic propagator for free 
Dirac fermions. 

The different forms of Xc and Xs lead to different re- 
sults for the full self-energy. We illustrate this for the 
T = case. For charge-mediated interaction, the soften- 
ing of the square-root singularity in Xs just implies that 



ImE(u,k) cx M 3 (9(e fc - uj) f 



C U 2 \u)\ 3 



(8) 



where c ~ vf/va, and f(x) subject to /(0) = 1 de- 
creases with increasing x and scales as 1/x at large x. 
This 1/x behavior implies that very near the mass shell, 
JmS(u, Ek) oc (ek — u>)9(ek — lj). The behavior of the 
spectral function for this case is shown in Fig. For 
experimental comparison, we added a constant impurity 
scattering 0. The spectral function is anisotropic with 
respect to Lu—tk, and has a kink at to—tk ~ u 2 u> 3 j ' E 2 F < 0, 
when the renormalization of the susceptibility becomes 
relevant. 





\ sy 


kink 














FIG. 2: The behavior of the MDC spectral function A(k,u) 
as a function of k at fixed uj for the interaction in the charge 
channel (a) and the spin channel (b). Observe that the kink 
in the spectral function is located on different sides of the 
mass shell in the two cases. 

For a spin-mediated interaction, the situation is dif- 
ferent. Due to the presence of the zero-sound pole in 
X s (0, g), ImS(u>, Eh) is non-zero at the mass shell, and 
only vanishes at ui — e k > ul\u>\ 3 /Ep > 0, when the 
Cherenkov-type absorption becomes impossible. The 
spectral function now has the form of Fig. [2). ft 
again has a kink (more precisely, a square-root non- 
analyticity), but now the kink is located at uj — ek > 0, 
i.e., on the other side of the mass shell. We see therefore 
that interactions in the charge or in the spin channel lead 
to qualitatively different results for the spectral function. 

A similar situation holds at a finite T. For breav- 
ity, we list the results for a charge-mediated interac- 
tion, and set k — kp. For the most interesting case 
T \u>\, we found that T 7 / 2 /y / |cJ| behavior holds down 
to \uj\ ~ oj = u 2 c T 3 /E 2 F . Below this scale, T 7 / 2 / y/\u~j~\ 
behavior crosses over to v/fjwf. The crossover behavior 
is somewhat involved, but the two limiting forms of ImE 
are captured by a simple extrapolation formula 

ImE( W> T)~u c |U(My *(x) = ^L. (9) 
Ef \ oj ) x+1 

We display this behavior in Fig|3] Note that at the max- 
imum, ImTi(u),T) scales as T 2 . This is qualitatively dif- 
ferent behavior from 7m£ oc T 3 , which would be the case 
if the mass shell singularity didn't exist. 
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FIG. 3: The schematic behavior of Im E(w, fcp,T) as a func- 
tion of a; at a given T. The dashed line is the result of the 
second-order in perturbation. 

The unusual behavior of the spectral function has 
a profound influence on the fermionic density of 
states N(oj) ~ f d 2 xImG(ui,x). In particular, 
N(0) oc T 7/3 . It, however, makes little differ- 
ence for the tunneling density of states N tun n{uj) oc 
(1/T) / duj'N{u')/ cosh 2 (to + u')/2T. The latter one is 
non-zero already in the non-interacting d— wave gas, 
where N(u>) oc \u>\. Then Nt unn (Q) oc T. The fermionic 
self-energy accounts for the corrections to the linear- 
T behavior. Typical frequencies and ek in the inte- 
gral for N tun n{0) are of the order of T, hence relevant 
JmS(u, k, T) are of the order of T 7 / 2 / VM ~ T 3 . As a 
result, N tunn /T = A + BT 2 is analytic in T 

To summarize, in this communication we re-considered 
the self-energy of a nodal fermion in a 2D d— wave 
superconductor. We found that the IrnE(u),€k,T) is 
actually a complex function of u),T and the quasi- 
particle energy tk, and is qualitatively different from 
ImT,(uj,T) oc mas(w 3 , T 3 ), suggested on general 
grounds. We found that the perturbative self-energy is 
non-analytic near the mass shell at T = 0. At a finite 
T, it diverges as T 7 / 2 / yf\uj\, when w tends to zero, and 
diverges even stronger near the mass shell. We demon- 
strated that the divergences in the self-energy for Dirac 
fermions originate from the same physics as the mass 
shell singularity in a 2D Landau Fermi liquid. We fur- 
ther demonstrated that the full summation of infinite di- 
agrammatic series eliminates the divergences and makes 
/mS(w,efc,T) finite, even at to = e^. However, the full 
S(ci),efc,T) at a given T is non-monotonic as a function 
of cj, and at a maximum is of order T 2 rather than T 3 . 

We also found that the full spectral function has a 
kink, whose location compared to the mass shell depends 
on whether the dominant fermion-fermion interaction is 
in the spin or in charge channel. For negative lo, relevant 
to ARPES experiments, kink is located at id — ej, > for 
spin-mediated interaction, and at id — tk < for charge- 
mediated interaction. We hope that careful ARPES stud- 



ies of the nodal spectral function in high T c materials 
below T c will be able to distinguish between charge and 
spin-mediated interactions. 

The doping dependence of the kink strength may also 
provide the information about the interaction. Singular 
self-energy along the nodal direction is determined by 
intra-nodal scattering, and the magnitude of the effect 
is determined by U(q — 0). This interaction is not en- 
hanced if lowering the doping drives the system closer 
to an instability at a finite q ( e.g., an antiferromagnetic 
instability). On the other hand, if the system evolution 
towards half-filling is governed by low-energy fluctuations 
at small q, one should expect a strong enhancement of 
the kink strength at smaller doping. 
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